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We investigate the dynamics of a wave packet in a parity-breaking one-dimensional periodic 
potential slowly varied in time and perturbed by a linear potential. Parity is broken by considering 
an asymmetric double well per unit cell. By comparing the prediction of the semiclassical dynamics 
with the full Schrodinger solution, we show that Bloch oscillations are strongly affected by anomalous 
velocity corrections related to Berry's phase. We characterize how these effects depend on the 
degree of parity breaking of the potential and on the modulation parameters. We also discuss 
how to measure the effects of the anomalous velocity in current experiments with non-interacting 
Bose-Einstein condensates in bichromatic optical lattices, under the effect of gravity. 

PACS numbers: 03.65.Sq,03.65.Vf,03.75.Lm 
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I. INTRODUCTION 

The semiclassical equations of motion play a funda- 
mental role in the transport theory for particles in crystal 
lattices, being them electrons in metals and semiconduc- 
tors [l|, or ultra-cold atoms in optical lattices Q. One 
of the most intriguing predictions, is the occurrence of 
Bloch oscillations when an additional constant force is 
present . Despite being hardly detectable for electrons 
in natural crystal lattices as predicted in Q , Bloch oscil- 
lations have been observed in a variety of systems, like 
electronic semiconductor superlattices |J| , optical wave, 
uide arrays Q , and ultracold atoms in optical lattices 

In the latter case, long-lived Bloch oscillations have 
been observed both with thermal clouds @ and with non- 
interacting degenerate fermions Q and bosons [|| loaded 
in sinusoidal periodic potentials under the effect of grav- 
ity, whose behavior is well described by the original semi- 
classical prediction ||. 

However, quantum mechanics has implications leading 
beyond the standard semiclassical approximation [l|, Q . 
In particular, for cyclic variations of the parameters, 
as is the case for Bloch oscillations in quasi-momentum 
space, phase dependent effects can be present, as shown 
in the seminal work by M. Berry Q. The importance 
of Berry's phase in the wave packet dynamics in peri- 
odic potentials was recognized by Chang and Niu [Ioj |. 
who derived a set of corrections to the semiclassical equa- 
tions of motions by means of an effective Lagrangian ap- 
proach [ll|. The whole set of corrections are now gen- 



erally called Berry corrections or Berry terms [12|. As 
a consequence, anomalous transport may occur in case 
of parity-breaking potentials, in two and three spatial 
dimensions. Anomalous effects can be present even in 
one-dimensional systems, provided that Bloch bands are 
also subjected to slowly-varying time-dependent pertur- 
bations, as discussed in [ll| (see also [l3l - [l6j ). 
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In this paper we investigate the effects of Berry correc- 
tions on the Bloch oscillations in one-dimensional (ID) 
time-dependent parity breaking potentials. In particu- 
lar, we consider the dynamics of a non interacting wave 
packet in a ID periodic potential, in the presence of an 
additional constant force, as can be obtained for instance 
with ultracold atoms in an optical lattice subjected to 
gravity 0, H|. Since in ID Berry terms may appear 
only when parity is broken and in the presence of time- 
depending bands [H, [l?} , we consider a potential com- 
posed by two shifted sinusoidal potentials with commen- 
surate wavelengths, whose amplitudes are modulated in 
time QJ]. 

We find that parity breaking affects significantly the 
band structure, already in the static case. Then, in 
the presence of time-modulations, Bloch oscillations 
may be dramatically modified owing to the presence of 
the anomalous velocity term in the semiclassical equa- 
tions. By direct comparison with the full solution of the 
Schrodinger equation, we show that these corrections of 
the semiclassical equations fully accounts for the center of 
mass dynamics of the wavepacket. We also characterize 
how the anomalous effects depend on the degree of par- 
ity breaking of the potential and on modulation parame- 
ters. Finally, we discuss how the effects of the anomalous 
velocity can be measured from time-of-flight measure- 
ments in current experiments with non-interacting Bose- 
Einstein condensates in optical lattices, in the presence 
of gravity. 

The paper is organized as follows. In Sect. HH we start 
by recalling the expressions for the semiclassical equa- 
tions of motions and the effect of Berry corrections on 
Bloch oscillations. Then, in Sect. IIII1 we discuss the 
model considered, and the static properties of the poten- 
tial in the presence of parity breaking. The emergence 
of an anomalous velocity term in the presence of time- 
dependent energy bands, its effect on Bloch oscillation, 
and the dependence on the system parameters, are dis- 
cussed in Sect. IIV1 Finally, in Sect. |V]we discuss the 
relevance of the present results for current experiments 
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with ultracold atoms in optical lattices. Conclusions and 
outlook are drawn in Sect. I VII 



II. ANOMALOUS VELOCITY 

Let us consider a quantum particle of mass m, repre- 
sented by a wavepacket ip(x), in a one dimensional peri- 
odic potential V(x) of period a, V(x + a) = V(x). Ow- 
ing to the Bloch theorem, the system eigenfunctions can 
be written as ip nk (x) — e u nk( x )i the u nk (xYs having 
the same periodicity of the potential. Consequently, the 
energies are restricted within energy bands e n (fc), that 
are periodic in the quasimomentum fc-space with period 
2ks = 2-ir/a, s„(k + k B ) = s n (k). 

At the semiclassical level, the particle dynamics can be 
described in terms of the wavepacket centers x c and k c 
in coordinate and quasi-momentum space, respectively. 
This holds even when the band structure is modulated in 
time, provided the variations are adiabatic. In particular, 
in the presence of an additional constant force F, and 
within the one-band approximation (we drop the index 
n), the evolution of the wavepacket centers is described 
by the semiclassical equations [l], [Io[ 
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where the Berry vector potential A k and scalar potential 
Xk in CQ are defined as [H [U] 
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A k (t) = «— <Ufc|-^litfc> , 
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(3) 

with the prefactor 2ir/a coming from the standard nor- 
malization of Uh^s in the unit cell, {u k \u k i) = (a/2iT)5kk' , 
at each time t [ljj . We recall that Zak-Berry phase 7 is 
defined as the integral over the first Brillouin zone of the 
vector potential A k [ItJ 
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According to Eq. ([2]), the quasimomentum center k c 
moves linearly across the first Brillouin zone, Tg = h/aF 
being the time needed to scan the full zone. The dynam- 
ics in real space is instead determined by the velocity 
term in the right side of in Eq. (JlJ, composed by the 
sum of the usual "normal" velocity 



v N (t) 
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and of the anomalous velocity 
dA k 



v A {t) 



dt 



dXk 
dk 



(5) 
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The latter represents the Berry correction to the standard 
equation of motions [l[ , and it appears as an electric field 




FIG. 1. Plot of the total potential [/(£) = V(£, 0) + 2£/t s at 
t = 0, for Vi = 10, V 2 = 5, 9 = 0.9tt. 



term in quasi-momentum space [2C 
locity can be cast into the form 



The anomalous ve- 



4-7T du k du k 
VA{t) = -- lla{ ^^k ) 



(7) 



that is manifestly vanishing for stationary potentials. In 
that case the wavepacket performs the usual Bloch oscil- 
lations, of period Tb and amplitude A/F, A being the 
energy bandwidth. Instead, in the case of time dependent 
parity-breaking potentials, the effects of the anomalous 
velocity may become relevant, as we discuss in the fol- 
lowing. 



III. THE MODEL 

In order to illustrate the effects of the anomalous veloc- 
ity in a specific case, we consider the following periodic 
potential 

Vx(t) cos 2 (qx) + V 2 {t) cos 2 (2qx + 9) (8) 



V(x,t) 

with period a = n/q. In the presence of an external force 
F, the hamiltonian takes the form 



h 2 d 2 



Fx 



(9) 



m being the particle mass. Here we will consider the 
specific case of gravity, F = mg. We also notice that the 
length q^ 1 , and the corresponding energy scale Er = 
h 2 q 2 /2m, represent two natural scales for the system, 
and it is therefore useful to rewrite Eqs. ([H])-© in a 
dimcnsionlcss form as 
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d 
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with 



V(& r) = Vi(t) cos 2 (£) + V 2 (r) cos 2 (2£ + 9) (11) 

where V; = Vi/Eji, £ = qx, r = Ept/%, and Tg = 
2itEn/mga is the dimensionless Bloch period, which cor- 
responds to the inverse of the normalized gravitational 
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FIG. 2. (Color online) Bandwidth A (a) and minimal 
bandgap E g (b), normalized to the gravitational energy per 
unit cell, as a function of 0, for Vi = V2 = 2 (red, continuous 
line) and Vi = 10, V2 = 5 (blue, dashed line). 



energy per unit cell, times 2n. The typical shape of the 
potential is shown in Fig. Q] Here we are mainly con- 
cerned with the single band approximation, and in the 
following we will refer to the (dimensional) bandwidth 
A of the lowest band as the bandwidth, and to the (di- 
mensional) minimal energy gap E g between the first and 
second band as the gap. 

As for the parameters, we choose typical values of the 
experiments with ultracold atoms in optical lattices Q. 
We fix tb/2tt = 10. Then, Vi is chosen in the range 
2 10 in order to have - at least in the case of a static 
primary lattice alone - clean Bloch oscillations (no drifts), 
large band gap (in order to conform to the single band 
approximation), and large enough oscillation amplitudes 
(i.e. large enough band widths A). The latter condition is 
obtained by considering Vi < 10. Conversely, the lower 
bound Vi = 2 guarantees that the ratio of the recoil 
energy to the gravitational energy is above the critical 
value for having no drift, tbVi/2tt « 20. In addition, 
for this value the minimal gap E g is already one order of 
magnitude greater than the gravitational energy per site 
mga (see Fig. and this accounts for the use of the 
single band approximation (actually, this corresponds to 
the adiabatic approximation, Tb > h/E g ). 

Let us now consider the effect of the secondary lat- 
tice. For V2 7^ 0, both the bandwidth A and the gap 
Eg strongly depend on 9, as shown in Fig. [2] [2l|. The 
band structure is periodic with period 71*, and is charac- 
terized by parity centers at integer multiples of 7r/2. In 
addition, A (E g ) also depend on the secondary lattice 
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FIG. 3. (Color online) Zak-Berry phase 7 as a function of 6, 
for Vi = V2 = 2 (red, continuous line) and Vi = 10, V2 = 5 
(blue, dashed line), as obtained from from Eq. Q. 



intensity, increasing (decreasing) monotonically as V2 is 
increased (in the range considered here) . In the following 
we consider the case V2 < Vi. 

For 9 7^ n7r/2 (n integer) parity is broken, and this is 
signaled by the appearance of a non vanishing Zak-Berry 
phase 7, defined in Eq. This quantity is proportional 
to the offset of the Wannier function centers UJ], u2l 



7 = — (w \x\wo) 



(12) 



(wn being the Wannier functions located at site j = 
[191 ]) modulo a vector of the direct lattice, coming from 
equivalent choices of the Bloch basis (with e 1 ^ being the 
gauge invariant quantity). The Zak-Berry phase 7 (cal- 
culated from Eq. (j4])) is shown in Fig. [3l as a function 
of 9, for two different sets of the potential intensities. In 
both cases, 7 grows monotonically reaching its maximum 
at 9 « 0.9-^0.957T, and then goes steeply to zero at 9 = n. 
The same behavior can be found in a wide range of values 
of Vi (i = 1,2). This suggests that the maximal effects 
of the anomalous terms in the dynamics should be ex- 
pected for these values of 9, as we will discuss in the next 
section. In addition, we find that the Zak-Berry phase is 
almost independent of Vi , and slightly decreases as V2 is 
increased. 

This analysis reveals that the control of the phase 9 is 
rather a crucial issue for commensurate bichromatic lat- 
tices, since the band properties have a strong dependence 
on it. This may become even more relevant in higher di- 
mensions, where parity breaking alone (no need for time- 
dependent bands) is sufficient for producing anomalous 
terms in the semiclassical equation of motions or other 
effects related to the Berry phase [l(| ■ 



IV. ANOMALOUS BLOCH OSCILLATIONS 

Let us now turn to the dynamical behavior of the sys- 
tem, by considering the evolution of an initial wavepacket 
under the effect of a constant force, and in the presence 
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of a modulation of the periodic potential. In particular, 
we consider a gaussian wavepacket modulated by the pe- 
riodic potential, obtained as the ground-state of the po- 
tential V(£, 0) plus an additional harmonic confinement 
V/io(£) = a £ 2 ( we cnose a = 3 • 10~ 5 , that corresponds to 
having about twelve occupied lattice sites) [23|. Then, 
the wavepacket evolution is investigated by solving the 
Schrodinger equation [24[ 



(13) 



with H(t) in Eq. ©. 

As anticipated, the combined effect of parity breaking 
and time modulation gives rise to an anomalous velocity 
term in the semiclassical equation for the center of mass 
dynamics in real space, see Eq. ([1]). The importance of 
this term can be analyzed by comparing the evolution of 
£ c (t) = (V'lClV , )r obtained from the full solution of the 
Schrodinger equation, with that predicted by the semi- 
classical equations ((I])-© 

In particular, let us separate the trajectory of the wave 
packet center into two parts 



(14) 



the first coming from the time integration of the "nor- 
mal" velocity in eq. (0 and the second from the "anoma- 
lous" velocity in eq. (JB]). Then, we are interested in 
studying the role of the contribution £f (r) on the full 
solution when the shape of the potential well is modu- 
lated in time, and to find regions in the parameter space 
where this term strongly characterizes the dynamics. We 
remind that Eq. © is not affected by anomalous terms, 
and the evolution of the wavepacket center in quasimo- 
mentum space is linear in time; in dimcnsionless units it 
reads 



k{r) = fc c (0) 



TB 



(15) 



Thus, we consider the following periodic modulations 
of the potential amplitudes 



Vi(r) = V 4 (l + ^sin 2 (fijr)) 



(16) 



where Ai and f2,; = r)in/TB are the modulation ampli- 
tudes and frequencies, respectively. 

Here we consider rational values of rjj, corresponding 
to periodic modulations with periodicity multiple of the 
Bloch period, in order to have a cyclic motion in param- 
eter (quasimomcntum) space, as required by the Berry's 
theory. The upper bound on rji, as well as that on the 
modulation amplitudes Ai, is given by the limit of valid- 
ity of the one band approximation, which requires that 
the band gap during the dynamics evaluated at k c (t), 
Eg(t), remains much greater than the modulation en- 
ergies. This condition turns out to be satisfied, in a 
large region of the parameter space, also for values of 
Ai,r)i 3> 1 [H| • Typical shapes of the potential during 
one modulation period are shown in Fig. 01 
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FIG. 4. Plot of the periodic potential V(£) at t/t b = 1/8 
(a), 1/4 (b), 3/8 (c), 1/2 (d), for Vi = 10, V 2 = 5, = 0.9tv, 
A!=A 2 = 5, 771 = 1, 772 = 2. 
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FIG. 5. (Color online) Evolution of the wave packet center 
for 741=^2 = 5, rji = 1, f]% = 2 (a), rji = 1, 772 = 3 (b), 
and 771 = 772 = 2 (c). The full semiclassical evolution (red, 
continuos line) perfectly fits the solution of the Schrodinger 
equation (red points). The terms corresponding to the normal 
dynamics, (black, dotted line), and anomalous dynamics, 
£^ (blue, dashed line), are also shown. 
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We start by describing some specific examples, shown 
in Fig. [5] Unless explicitly stated, all the results pre- 
sented in this section are obtained for the following set 
parameters characterizing the "static" properties of the 
potential: Vi = 10, V 2 = 5,0 = 0.9tt. Then, we will 
discuss the dependence of the results on the choice of 
parameters, as well as their generality. 

In Fig. [5]we plot the evolution of the wave packet cen- 
ter £, c (t), along with the separate contributions coming 
from the normal and anomalous terms, for Ai = Ai = 5, 
Vi = 1) f\i = 2 (a), i]i = 1, r] 2 = 3 (b), and 771 = 772 = 2 
(c). The full semiclassical evolution of Eq. ([1]) (red, 
continuos line) perfectly agrees with the solution of the 
Schrodinger equation (red points). For r/i ^ 772, the nor- 
mal term £^(t) (black, dotted line) shows peaks only at 
t = and t = tb whereas the anomalous term £;f (r) 
(blue, dashed line) shows one (a) or two (b) additional 
central peaks. In general, we find that the number of cen- 
tral peaks (excluding those at r = and r = tb) equals 
i]2 — 1. This holds also for higher integer values of 772, 
provided that the single band approximation is satisfied. 
The dependence on 771 at fixed 772 is more weak, affecting 
the trajectory without basically changing its structure. 
In particular, when 771 = 772 as in Fig. [5]:, the modulation 
does not change the shape of the potential (if A\ = A2). 
In this case we find an additional central peak also in the 
normal term, in the opposite direction with respect to 
the anomalous component. 

We notice that in all these examples the anomalous 
term, Cc ( t )j dominates over the normal one, £^(t). In 
general, the relative weight between the two depends on 
the amplitudes A i: as will be discussed later on. 

We have also analyzed cases in which 771 , 772 < 1 , find- 
ing that the qualitative picture does not change. The 
significant difference is that the periodicity of the motion 
amounts to several Bloch periods, as expected. In ad- 
dition, even larger values of Ai can be choosen without 
breaking the single band approximation. 

An interesting issue concerns the dependence of these 
dynamical results on 9, which can be considered as a de- 
pendence on the degree of parity breaking. In Fig. [5] 
we plot the evolution of the wave packet center £ c (r) by 
varying 9, for A\ = A2 = 5, 771 = 1 and 772 = 2. In panel 
(a) we report the full solution £ c (t), while in (b) just the 
contribution of the anomalous term £f (r). The values 
of 9 correspond to a growing Zak-Berry phase (see Fig. 
[3]) from 7 = (no parity breaking) at 9 = 0.5n, towards 
its maximum value around 9 = 0.9ir. It is evident that, 
as anticipated in the previous section, the amplitude of 
the anomalous oscillations grows together with 9 in this 
interval, showing the relevance on the degree of parity 
breaking on the resulting dynamics. It is worth men- 
tioning that in the interval of negative Zak-Berry phases, 
that is for 9 G (0,7r/2) the behaviour is almost specular. 
Indeed, the two cases are connected by a mirror transfor- 
mation 6 — > 7T — 6; though this is not an exact symmetry 
due to the presence of gravity, the qualitative picture 
does not change in the two cases. 
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FIG. 6. (Color online) (a) Evolution of the wave packet center 
£ c (t) for A\ — A2 = 5, 771 = 1,772 = 2 and different values of 
8. (b) The contribution of the anomalous term £f (t). 



Let us now consider the dependence on the modulation 
amplitudes Ai. Since the anomalous velocity contains 
time derivatives of the periodic Bloch wave functions, 
we expect the anomalous effects to increase as Ai is in- 
creased. This is indeed the shown in Fig. [7]where 
we plot the full solution £ c (t) (a) and just the anoma- 
lous term £f (t) (b), for different values of A = A\ = A2, 
and 771 = 1, 772 = 2. This figure shows that, for A in 
the range [0, 10], the amplitude of anomalous oscillations 
grows monotonically as expected (we find that in gen- 
eral higher values of A are beyond the validity of the 
single band approximation). In addition, we find that 
while the amplitudes of normal Bloch oscillations (for 
^4 = 0) can be of the order of some lattice site, those of 
anomalous oscillations are rather smaller. In particular, 
anomalous oscillations become predominant only when 
the bandwidth is small (that is, for rather large Vi(t)). 

Finally, we have also studied the case A\ ^ A2, find- 
ing a more pronounced dependence of the results on A2 
rather than on A\, according to what previously found 
for 771 and 772. This is not surprising, since parity is bro- 
ken due to the presence of secondary lattice, and conse- 
quently its parameters mainly influence the anomalous 
dynamics. 



V. MEASURING THE ANOMALOUS 
VELOCITY 

Let us now discuss the experimental relevance of these 
results, focusing on current experiments with ultracold 
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FIG. 7. (Color online) (a) Evolution of the wave packet center 
£ c (t) for r)i = 1,772 =2 and different values of ^4 = Ai = Ai- 
(b) The contribution of the anomalous term £^(t). 
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FIG. 8. (Color online) Plot of the momentum distribution 
|*(p)| 2 at t/t b = (a), 1/4 (b), 1/2 (c), 3/4 (d), for the case 
shown in Fig. [5ji. The horizontal axis on the top represents 
the distance after the time-of-flight (see text). 



atoms in optical lattices. As we have previously dis- 
cussed, the amplitude of anomalous Bloch oscillations 
in coordinate space is rather small (only a fraction of 
a lattice size), and therefore a direct detection by in- 
situ imaging is not feasible, due to resolution limita- 
tions. Usually, Bloch oscillations are detected via time- 
of-flight measurements, probing momentum space. Typ- 
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FIG. 9. (Color online) The wavepacket velocity for the case 
shown in Fig. [SJi. (Ai = A2 — 5, 771 = 1, 772 = 2). The semi- 
classical velocity ^ c (red, continuos line) fits that extracted 
from the momentum distribution of the full Schrodinger so- 
lution (red points). The normal velocity ^ is also shown 
(black, dashed line). The right vertical scale corresponds to 
the center (x)t elp of the density distribution after a time of 



flight t e 



50 ms for a non-interacting condensate of K. 



ical shapes of the momentum distribution ^(p)! 2 are 
shown in Fig. [5] (for the wavepakct in Fig. [S^,) . It 
is characterized by sharp peaks at p(r) = ±2n + k c (r) 
(n = 0, ±1, . . .) [H, H3' with k c (r) being the wavepaket 
quasimomcntum given by Eq. (|15|) . This equation is 
not affected by anomalous terms, and therefore it is not 
possible to measure the effect of the Berry terms sim- 
ply from the shift of the peaks, as they move linearly 
in time as for normal Bloch oscillations. Nevertheless, 
their relative population depends on the modulation, and 
the wavepacket velocity £ c (r) can be extracted from the 
momentum distribution, according to Ehrenfest theorem 
d{x) /dt = (p) /to, that in dimensionles units reads 



Ur) = 5<P>(r). 



(17) 



In the experiments, the momentum distribution is in- 
ferred from the density distribution after the time of 
flight , owing to the relation 



\4>oo(x,t)\ 



lim \ip(x,t;t e )\ 2 oc \t>(mx/t e ,t)\ 2 (18) 

t e — >oo 



where t e is the free expansion time, and t the evolution 
time in the trap. Therefore, the wavepacket velocity £ c (t) 
after an oscillation time r = Ent/h can be measured from 
the center (x) tc = {ip(t;t e )\x\ip(t;t e )) of the measured 
density distribution after a time of flight t e , as 



; n 2 m / \ 



(19) 



In Fig. [9] we show the comparison of the semiclassi- 
cal velocity £ c (t) with that extracted from the momen- 
tum distribution, 0.5 (p)(r), and the corresponding cen- 
ter (x)t e of the density distribution after a time of flight 
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FIG. 10. (Color online) Plot of maximum of the modulus of 
the normal and anomalous velocities, |^ | (empty symbols) 
and (filled symbols) respectively, as a function of 9, for 
T)i = 1, r]2 = 2 (squares) and r?i = 1, r\% = 3 (circles). In all 
cases Ai = A2 — 5. 

i e = 50 ms, for a non-interacting 39 K condensate In 
this case, (x) te turns out to be of the order of some tens 
of microns, and its detection is experimentally feasible. 
In addition, the effects of the anomalous velocity term 
are pronounced in a large time range where the normal 
velocity is almost vanishing, and this should help their 
visibility. 

As discussed before, these effects can be maximized by 
chosing suitable values of the angle 6. In Fig. [TU]we plot 
the maximum of the anomalous and normal velocities 
and ^ , respectively) as a function of 9, for 0.57T < 
9 < 7r and two different sets of rj^s. In particular, it is 
convenient to choose 9 ps 0.97T, where the magnitude of 
£^ is close to its maximun and ^ almost vanishing. 



VI. CONCLUSIONS 

We have investigated the occurrence of anomalous 
Bloch oscillations in the dynamics of a wave packet in a 
parity-breaking one-dimensional periodic potential that 
is modulated in time, in the presence of an additional 
weak linear potential. We have considered a periodic po- 
tential composed by the sum of two sinusoidal potentials 
of commensurate wavelengths, shifted by an angle 9 that 
characterizes the parity breaking. We have found that 
the anomalous velocity correction of the semiclassical 
equations - related to Berry's phase - profoundly affects 
the wave packet dynamics, in agreement with the full so- 
lution of the Schrodinger equation. We have investigated 
how these effects depend on 9, finding that the maximal 
effects are for 9 « 0.97T both for static and dynamical 
properties. These results are an explicit demonstration 
of the importance of Berry's corrections beyond the usual 
semiclassical approximation. In particular, these effects 
may become relevant for current experiments with ul- 
tracold atoms in optical lattices, when also time mod- 
ulations and parity breaking are considered. As an ex- 
ample, we have shown that the effects of the anomalous 
velocity can be measured from the atomic distribution 
after the time-of- flight, in current experiments with non- 
interacting Bose-Einstein condensates in bichromatic op- 
tical lattices. 
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